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ABSTRACT 

We investigate the impact of non-Gaussian scatter in the cluster mass-observable scaling relation 
on the mass and redshift distribution of clusters detected by wide area surveys. We parameterize 
non-Gaussian scatter by incorporating the third and forth moments (skewness and kurtosis) into 
the distribution P {Mohs\M). We demonstrate that for low scatter mass proxies the higher order 
moments do not significantly affect the observed cluster mass and redshift distributions. However, 
for high scatter mass indicators it is necessary for the survey limiting mass threshold to be less than 
10^'*/i~^Mq to prevent the skewness from having a significant impact on the observed number counts, 
particularly at high redshift. We also show that an unknown level of non-Gaussianity in the scatter 
is equivalent to an additional uncertainty on the variance in P [Mahal M) ^^'^ thus may limit the 
constraints that can be placed on the dark energy equation of state parameter w. Furthermore, 
positive skewness flattens the mass function at the high mass end, and so one must also account 
for skewness in P(A/obs|-^^) when using the shape of the mass function to constrain cluster scaling- 
relations. 

Subject headings: cosmology: dark matter — galaxies: clusters: general — intergalactic medium 



1. INTRODUCTION 

The evolution of the number den sity of galaxy clusters 
is a sensitive c osmological probe (jBahcall fc FanI [79981 
lEke et aLlll998f ). As an indicator of the expansion rate 
as a function of time, the galaxy cluster number den- 
sity is sensitive to t he dark energy eq uation of state 
(jHaiman et all 1200X1: \We\leT et al.ll2001h . This provides 
a growth-based dark energy test, an important comple- 
ment to the distance-based tests that have provided the 
most compelling evidence for dark energy t o this point 
(jPerlmutter et al.lll999HSchmidt et al.lll998f) . 

Galaxy clusters can be selected by many diverse meth- 
ods, including (but not limited to) optical richness. X-ray 
thermal bremsstrahlung flux, weak lensing shear and the 
Sunyacv-Zcl'dovich (SZ) effect. The key challenge for 
using galaxy clusters as precise cosmological probes is 
in understanding how to relate observables to a quan- 
tity that can be well predicted by theory, for example, 
mass. The ultimate goal is to produce theoretical pre- 
dictions of the distributions of observables as a function 
of redshift and cosmological parameters. Short of this, 
one approach is to theoretically model the evolution of 
number density as a function of mass, and then estimate 
the mapping between observables and mass in order to 
predict the observed evolution. This mapping can either 
be estimated from theoretical considerations or be de- 
termined directly fr om the data, assuming some regular- 
ity in the mapping (fTyTaiumdar fc Mohij 120031: IH^|2003I : 
ILima fc Hull20M . ^0m 

It is important to understand the statistics of the rele- 
vant mass-observable scaling relation. Although clusters 
may follow a mean relation, individual clusters will de- 
viate from it. If the level of scatter around the mean 
relation is not small, the shape and amplitude of the ob- 
served mass function can change significantly. At cluster 
scales the mass- function is a steeply declining function of 
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mass. Therefore a larger number of low mass clusters will 
scatter over the detection threshold of the survey than 
those higher mass clusters that scatter below it. The 
net increase in the total number of clusters in the sam- 
ple is thus dependent on the slope of the mass function 
at the threshold mass, and the magnitude of the scat- 
ter around the mass-observable scaling relation. If the 
latter is well constrained then the measurement of the 
cluster mass-function is actually improved statistically 
due to the reduction in shot noise. However, in practice 
it is difficult to precisely measure the scatter and the- 
oretical estimates can vary substantially. Furthermore, 
as we shall demonstrate, for large intrinsic scatter, the 
magnitude of the higher order moments (skewness and 
kurtosis) can become significant. 

Previous work forecasting the constraints on cosmolog- 
ical parameters that can be achieved by cluster surveys 
have assumed the scatter in scaling relations to be log- 
normally distribution around the mean relation (that is, 
normally distributed in the logarithm of the mass). How- 
ever, cosmological simulations and observations of large 
samples of clusters have demonstrated deviations from 
lognormal behaviour. In general, the cause of these de- 
viations can be seperated into two catagories; dynamical 
state and projection effects. 

The former refers to the impact of a sub-population 
of clusters that are systematically offset fro m the mean 
scaling rela tion due to the i r dyn amical state. I Shaw et al.l 
(|2006) and lEvrard et all (|200^ demonstrated using N- 
body simulations that dark matter halos undergoing a 
major merger have a systematically higher dark matter 
velocity dispersion (udm) t han their more relax ed coun- 
terparts (of the same mass). lEvrard et al.l (|2008[ ) showed 
that the small fraction of interacting halos cause posi- 
tive skewness in the distrib ution of the r e sidual s around 
the mass — (Tdm relation. IStanek et al.l ()2009[ ) analyse 
the covariance of bulk cluster properties using a large 
sample of halos extracted from hydrodynamical Simula- 
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tions. They find small deviations from Gaussian scatter 
for some properties, most nota bly in the mass-we ighted 
temperature ~ mass relation. iPratt et al.l (|2009l ) anal- 
ysed X-ray luminosity scaling relations using a sample of 
31 nearby clusters observed by the XMM-Newton X-ray 
observatory. They find that the scatter in the L^ — Yx re- 
lation (where Yx was assumed to be a robust, low scatter 
proxy for cluster mass) was significantly non-Gaussian 
due to the systematically above-average luminosity of 
cooling core clusters in their sample. They also note 
that the scatter becomes more Gaussian when the cen- 
tral regions of the clusters are excluded in the luminosity 
measurements . 

The second cause of non-Gaussianity in observed scal- 
ing relations is confusion in cluster selection due to pro- 
jection effects. Several studies have demonstrated that 
lower mass, unresolved clusters, as well as gas outwith 
cluster regions can contribute significantly to the mea- 
sured integrated SZ flux (Y), causing additional scat- 
ter in the Y — M relation and introducing a tail to- 
wards high flux in t h e distribution of Y a t constant mass 
(I White et all l2002t iHolder et al] l2007t iHallman et all 
|2007| ). A similar e ffect is found for optically-selected 
clusters. ICohn et~a l. (2007) demonstrated using mock 
galaxy catalogues that a non-negligable fraction of clus- 
ters (10% at z = 0.4 to 22% at z = 1) identifled using the 
red-sequence are 'blends' - cluster candidates in which a 
large number of different halos have contributed galax- 
ies. Blends thus cause a tail towards high richness in the 
distribution of optical richness at fixed mass. 

In this work, we relax the assumption of lognormal 
scatter and investigate the effect of non-Gaussian scat- 
ter around the mean mass-observable scaling relation on 
the observed mass and redshift distribution of clusters. 
Specifically, we quantify the impact of non-zero skew- 
ness and kurtosis - the third and fourth standardized 
moments of a generalized probability distribution - as 
perturbations to the purely Gaussian case. 

Throughout this paper we make a distinction between 
true mass M and observed mass Mobs- The former is 
the actual cluster mass, as defined by a spherical over- 
density, A, and meas ured in cosmological simulations of 
structure formation (iJenkins et al.1 120011 : I Warren et al.l 
120061 iTinker et al.ll2008l ). M„h„ is the cluster mass that 
would be inferred from observations through application 
of a scaling relation (e.g. Y — M, — M, etc), or via 
self-calibrati on of the mass function ()Lima IFHuI [20051 : 
ICunhal [20091) . 

2. SCATTER IN THE MASS-OBSERVABLE SCALING 
RELATION 

The predicted redshift distribution of clusters observed 
by a given survey is given by 



dN _ dV dM dn 

d7 ~ dzdn Jo ~M"dln(M) 



f{M,z), (1) 



where M is cluster mass and dn/d\nM is the mean co- 
moving number density of clusters (the mass function). 
The function /(Miim, z) represents the survey selection 
function which accounts for the limiting mass of the sur- 
vey, Afiim defined by some threshold in in the mass ob- 
servable (e.g. optical richness. X-ray luminosity, inte- 
grated SZ flux), and the statistics of the mapping be- 
tween Mobs and M. In the limit of perfect (zero scatter) 



mass measurements, this is simply a step function at the 
limiting mass of the survey. 

If one assumes lognormal scatter around the mean scal- 
ing relation (Gaussian scatter in InAf) then the proba- 
bility P (Afobs|M) of observing the mass Afobs given the 
"true" underlying mass M, is 
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This parameterization allows for redshift dependent scat- 
ter CTinM and bias Afbias in the rnass-ob servable scaling 
relation ([Lima fc Hull2005HCunhall2009( l. For simpHcity, 
we henceforth ignore the bias term and concentrate on 
the impact of scatter alone. The distribution of observed 
cluster masses is just a convolution of the true mass func- 
tion with P(Afobs|Af), 



dn 



dlnA/obs 
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d\nM 



P(Afobs|Af)dlnA/ . (4) 



Plugging in Eqn. [2] and assuming an intrinsic power law 
distribution in mass, dn/d\n{M) oc M~", it is straight- 
forward to show that the observed mass distribution is 
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where o denotes the true mass distribution evaluated at 
Afobs by applying the mean M — Mobs scaling relation. 

The extent of the deviation of the observed mass func- 
tion from the true mass function is clearly controlled by 
Q!crin M , the product of the standard deviation of the dis- 
tribution P(Mobs|M) with the slope of the mass func- 
tion, a. Evidently, a constant acinM results in an ob- 
served mass function that has a constant and positive 
vertical offset from the true mass function. For more re- 
alisti c mass functions ([Jenkins et al.l 120011 : [Tinker et al.l 
[2008f ) the slope a is approximately 1 at the group mass 
scale and increases with increasing mass and redshift, 
exponentially so at very high mass and redshift. The 
impact of scatter on the observed mass function is thus 
significantly greater at high masses/redshifts (see Section 

Above a limiting threshold in Afobs, scatter in a clus- 
ter scaling relation causes a net increase in the number 
of detected clusters. An unknown amount of scatter in 
the scaling relation thus degrades the cosmological con- 
straints that can be obtained from cluster nu mber counts 
as on e must also marginalise over crinM- [Lima fc Hu[ 
([2005( ) demonstrated that for a survey with a fixed limit- 
ing mass of 10^^-^/i~^Mq and (rf^^ = 0.25^, a la uncer- 
tainty of 0.25^ on af^ would produce a 10% uncertainty 
in the number counts a,i z = 0.5, a 20% uncertainty at 
z = 1, and a 50% uncertainty at z = 2. 

2.1. Non-Gaussian Scatter 

We now determine the impact on the observed 
mass function of non-Gaussian scatter in the mass- 
observable relation. We proceed by using the Edge- 
worth series to approximate a non-Gaussian P (Mobs|M) 
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(e.g. Bernardeau fc Kofmanlll995tlBlinnikov fc Moessneil 



I1998D . thus taking it to be a perturbation to the Gaus- 
sian case. It is parameterised by the third and fourth mo- 
ments, the skewness (7) and kurtosis (k), of the probabil- 
ity distribution P {Mohs\M). The Edgeworth expansion 
is particularly useful for convolutions when expressed as 
a series of derivatives of Gaussians, 

, , 7 d^G K d'^G 72 d^G , , 
P(M.gM).G(x)-|^-H-^ + 2^^ (6) 

where the skewness, 7 is defined as 

^3\ 

7 = 3 (7) 

and the kurtosis, k as 



((Afobs-Af)') 



(A^ob. - A//)' 



(8) 



and G{x) is a Gaussian distribution in x (equal to 
P(A^obs|M) inEq.©. 

By plugging Equation [6] into Equation [4] and doing 
some integration by parts, the observed mass function 
can be calculated for any true mass function, 
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Assuming again an intrinsic power-law distribution of 
mass with slope a, dn/dln A/truc — Af~", it is straight- 
forward to show that the observed mass distribution is 



now 
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where o denotes the true mass distribution evaluated at 

Afobs. 

The relevant parameter here is again clearly aainM- 
If this parameter is large, there are two important ef- 
fects: the number of objects at any given mass scale is 
increased substantially (as for the Gaussian case), and 
the higher order moments of the distribution, 7 and k, 
become important. The transition at which the latter oc- 
curs is clearly when aainM becomes greater than unity. 
For example, if acrinM = 1, then 7 = 1 produces an 17% 
increase on the number counts compared to the purely 
Gaussian case at any given mass scale. Assuming k — 1 
provides an additional 4% correction. Note that Equa- 
tion [10] demonstrates that both (positive) 7 and k cause 
an up-scattering of clusters. For skewness this is due to 
the tail towards large mass that increases the probability 
of a cluster of having Afobs >> A/. For kurtosis, the up- 
scattering is due to the wider wings of the distribution. 

3. IMPLICATIONS OF NON-GAUSSIAN SCATTER FOR SZ 
AND OPTICAL CLUSTER SURVEYS 
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Fig. 1. — Contours of constant slope, a of the halo mass func- 
tion, where dn/dln{M) oc M~°', as a function of mass and redshift 
for our fiducial cosmology. The solid black contour emphasizes 
the critical curve Mc{z) (see text) for a DES-like optical survey 
(""inM = 0-5)i the black dashed contour shows Mc(z) for an SPT- 
like SZ survey {cri^m = 0.25). The black solid and dashed vertical 
lines show the assumed limiting masses for the DES and SPT sur- 
veys, respectively. 



We now evaluate the impact of non-Gaussian scat- 
ter in the mass-observable scaling relation on the pre- 
dicted mass and redshift distribution of clusters observed 
by a South Pole Telescope-like SZ survey and a Dark 
Energy Survey-like optical survey. For the SZ survey 
we assume that th e scatter in the m ass-SZ fiux rela- 
tion is (Tsz = 0.25 (|Shaw et al.l |2008[ ) and the hmiting 
mass is 3 X 10^''/i~^Mo. For the optical survey we as- 
sume the scatter in the optical richness - mass relation 
is (Topt = 0.5, in keeping with the results of Rykoff et "aD 
(2 008a.h ); Becker et al. (2007); Rozo et al. (2009, 200^), 
and a constant limiting mass of 5 x lO^'^/i~^M0. Note 
that iBecker et al.l (|2007f) found that the scatter in the 
richness-mass relation is 0.5-0.75, and we have taken the 
value at the lower end of this range. 

We have demonstrated that the impact of the higher 
order moments 7 and k, become significant when fji^na, 
the product of the slope of the mass function and the 
standard deviation of P {Mohs\M), becomes greater than 
unity. Under the simplifying assumption that crinM is 
independent of mass and redshift for a given observable, 
the higher order moments thus become relevant when the 
effective slope of the mass function exceeds the threshold 
value, Uc > l/cinM- This critical slope is reached above 
a redshift (and cosmology) dependent mass threshold 
Af,(z). For SPT, ac,sz = 1/0.25 = 4, for DES ac,opt = 2. 

In Figure [T] we plot the slope of the mass function 
as a function of mass and redshift . For t his, we adopt 
the mass function of iTinker enil ()2008[ ). and assum- 
ing dn/d\n{M) oc M~" at any given mass and redshift, 
calculate a over a wide range of mass and redshift. The 
mass function is defined in terms of AI20Q, where the sub- 
scipt denotes the mass enclosed in the region of spherical 
overdensity of 200 times the mean density of the Universe 
(at the relevant redshift). For our fiducial cosmology we 
assume parameters consisten t with the WMAP 5-year 
results (Qm = 0.27, erg = O.S.'Dunkl ev et al.ll2009[ ). The 
grey lines denote contours of constant a, from 10 (top 
right corner) to 1 (bottom left corner). The thick black 
solid and dashed contours denote a = 2 and 4, the critical 
slopes ttc for DES and SPT, respectively. These contours 
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Fig. 2. — (Upper) The observed redshift distribution of clusters, 
dN/ dz above a limiting observed mass threshold of 5 X IO^^/i-IMq 
and 3 X lO^'*/i~'^M0. (Lower) The observed mass distribution of 
clusters M'^^^/ pra(dn/ dMahs) at z = 0.25 and z = 1. In each 
plot, the line colors represent the value of crinM^ (blue), 0.25 
(black) and 0.5 (red). The solid lines represent 7 = k = 0, the 
dashed 7 = 0.5, k = 0, dot-dashed re = 0.5 (7 = 0) and the dotted 
line 7 = 1, K = 1 (in our fiducial cosmology). The blue dotted 
line shows the results for a WMAP3 cosmology, where erg = 0.77, 
f^M = 0.24, Ho = 73km/s/Mpc, and we take (TinM = 0. 

thus give Mc{z) for each survey. The vertical soHd and 
dashed Hues represent the hmiting mass for each survey. 

The Figure demonstrates that for the SZ survey, aa\n m 
becomes greater than unity above 8 x lO^''/i~^M0 at z= 
0.5, and above 4 x lQ^^h~^MQ at z = 1. Given the 
scarcity of objects above these thresholds, it is clear that 
7 and K, will not significantly impact the SPT cluster 
number counts. 

For the optical survey, Mc{z) is 3 x lO^^/i^^M© at z = 
0.2, 1.7 X IQi^/i-iMo at z = 0.5 and 5.25 x IO^^/i-iMq 
at z = 1.2. As Mc{z) remains greater than the expected 
limiting mass for DES, 7 and k will not strongly effect the 
total number of detected clusters. However, a sizeable 
number of clusters will be detected in mass bins that 
exceed Mc{z). Non-Gaussian scatter in the mass-richness 
relation may therefore cause a detectable redistribution 
of clusters in mass bins greater than Mc(z), 

In Figure [2] we demonstrate the impact of skewness 
and kurtosis in P (Mobsl-^^) on the observed redshift and 
mass distribution of clusters. In both panels, the colors 
of each line represent the value of ainM; (blue), 0.25 
(black) and 0.5 (red). The solid lines represent (in our 



fiducial WMAP5 cosmology) 7 = k = 0, dashed 7 = 
0.5, K = 0, dot-dashed k = 0.5 (7 = 0) and dotted 7 = 1, 
K = 1 (as a more extreme scenario). The blue dotted line 
shows the results for a WMAP 3-year cosmology, where 
0-8 = 0.77, VLm_= 0.2 4, Hq = 73km/s/Mpc, and a\nM = 
(|Spergel et al.l[2007i) . 

In the upper panel we plot the redshift distribution of 
clusters dN/dz using our two A/obs limits; 5 x 10^'^/i~^Mq 
for DES (upper lines) and 3 x lO^^h'^MQ for SPT (lower 
lines). For the purposes of comparison, we have plot- 
ted the results for CTinM = 0.25 and 0.5 for both hmit- 
ing masses. As predicted by Equation [TOl the impact of 
skewness and kurtosis on the redshift distribution is de- 
pendent on the value of din m and the local slope of the 
(true) mass function at the limiting mass of the survey. 
Thus at higher redshifts the impact of 7 and k becomes 
greater. The flattening of dN/dz due to the higher or- 
der momements is more pronounced than that due to an 
increase in erg or cinM- 

For a threshold of Mobs.iim = 5 x 10^^/i"^Mq and 
cinM = 0.5 a skewness of 7 (k) = 0.5 increases the 
mean number of clusters by 10% (2%) at 2 = 0.75. The 
more extreme case of 7 = k = 1 produces a 24% in- 
crease. For the same scatter but with a limiting Mobs 
of 3 x 10^'^h~^MQ, 7 (k) = 0.5 increases the number of 
clusters at z = 0.75 by 42% (14%), and 7 = k = 1 by 
218%. It is clear that it is the low mass threshold of 
the optical survey that prevents the higher order mo- 
ments from strongly affecting the number of detected 
clusters in each redshift bin. The increase in dN/dz (at 
z = 0.75) for 7 — 0.5 is roughly equivalent to assuming 
lognormal scatter but increasing a\nM by 22% to 0.61 
(for Mobs,iim = 3 x IO^-^/i-^Mq) and by 14% to 0.57 
(Mobs,iim = 5 x IO^^H^^Mq). For our extreme case of 
7 = K = 1, the results are equivalent to increasing fiinM 
to 0.73 and 0.65 at each mass threshold, respectively. 

For the SZ survey (cri„M = 0.25, Mobs.iim = 3 x 
10^'^h~^MQ), increasing 7 from to 0.5 provides a 8% 
increase in the number counts at z — 0.75. This is equiv- 
alent to an increase in ain m for purely lognormal scatter 
of 12% (to 0.28). The same change in kurtosis provides 
less than a one percent change in the number counts. 
The 7 = K = 1 case provides a 20% increase in the num- 
ber counts, equivalent to an increase of cinM by 28% to 
0.32. 

In the lower panel of Figure [5] we plot the mass distri- 
bution of detected clusters, M'^^^/ p^{dn/ dMohs) (scaled 
to reduce the dynamic range, where pm is the mean mat- 
ter density of the Universe at each redshift) at z — 0.25 
and 1.0. The line types and colors represent the same 
values of crinM, 7, and k as the upper plot. For inte- 
grated SZ flux (for which (TinM = 0.25), the higher-order 
moments do not have a significant effect on the shape 
of the mass function. The impact of 7 and k is more 
evident for high-scatter mass indicators (red lines), flat- 
tening the mass distribution above 10^'*'^/i~^Mq at z — 
0.25, and 10^''/i~^Mq at z = 1.0. These masses roughly 
correspond to the predicted values of Mc(z) at these red- 
shifts in Figure [1] (a = 2 contour) . 

Overall, Figure[2]demonstrates two points. A skewness 
of 7 = 0.5 has the same effect on dN/dz (at z = 0.75) 
as increasing the variance (Ji^m by 0.075 for DES and by 
0.016 for SPT. For the significantly more non-Gaussian 
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distribution parameterised by 7 = k = 1, the equiva- 
lent change in erf is 0.17 for DES and 0.04 for SPT. 
iLima &: Hul (|2005D show that a prior of a-{<7f^M) = ^-Ol is 
necessary to ensure that uncertainty on the level of scat- 
ter does not degrade the constraints that can be placed 
on the dark energy equation of state parameter w (see 
their Figure 6). Therefore, assuming 7 does not greatly 
exceed 0.5, non-Gaussian scatter in the mass-SZ flux scal- 
ing relation should not significantly degrade the cosmo- 
logical constraints that can be achieved by SPT. How- 
ever, for surveys utilizing a high-scatter mass proxy such 
as optical richness, an unknown level of non-Gaussianity 
in P {Mohs\M) may limit the constraints that can be 
placed on the dark energy equation of state via cluster 
number counts. 

Secondly, if one aims to use the shape of the mass 
function to self-calibrate a survey to obtain information 
on the slope, normalization and variance of the mass- 
observable scaling relation, then care must be taken to 
ensure that the parameters are not biased due to the 
increasing impact of the higher order moments on the 
shape of the mass function towards high masses. One 
means of ensuring this is to include only a single mass 
bin for Mobs > Mdz). 

4. CONCLUSION 

We have investigated the impact of non-Gaussian scat- 
ter in the mass-observable scaling relation on the cluster 
mass function by incorporating the third and forth mo- 
ments (skewness and kurtosis) into the probability dis- 
tribution P (Mobs I A^) via the Edgeworth expansion. We 
have demonstrated that if aainM - the product of the 
standard deviation in P {Mohs\M) and the slope of the 
mass function at the limiting mass of a survey - is greater 
than unity, then positive skewness and kurtosis will in- 
crease the number of clusters detected and flatten the 
high-mass end of the observed mass function. 

For low scatter mass proxies like integrated SZ flux, 
higher order moments do not significantly affect the mass 
and redshift distribution of clusters. However, for sur- 
veys utilizing a high scatter mass proxy such as optical 



richness, the limiting mass threshold must be less than 
10^'*/i~^Mq to ensure that the skewness does not signifi- 
cantly effect dN/dz, especially at high redshift. We have 
also found that an unknown level of non-Gaussian scat- 
ter is roughly equivalent to an additional uncertainty on 
the variance ctj^^m ^ (-^^obsl^f) and thus may limit the 
contraints that can be placed on the dark energy equa- 
tion of state parameter w, especially for surveys that use 
a high scatter mass proxy. Furthermore, if one wishes to 
use the shape of the mass function in each redshift bin to 
self-calibrate for cluster scaling relation parameters then 
it will be necessary to account for non-Gaussian scatter 
on the shape of the mass function by marginalising over 
the skewness and kurtosis parameters, 7 and k in addi- 
tion to the the variance ^ , slope and normalisation in 
each redshift bin. 

We note that the values of 7 and k used in the ex- 
amples given in this work were chosen arbitrarily (and 
to be small enough to ensure that the Edgeworth ex- 
pansion remains an appropriate approximation to a non- 
Gaussian P{Mohs\M)). Large- volume simulations and 
mock-galaxy catalogues containing several thousands of 
clusters will be necessary to obtain better motivated pre- 
dictions of 7 and K for optical, SZ (and X-ray) surveys. 
Measuring 7 and k would require a large observational 
sample with precisely measured masses (for example, us- 
ing X-ray spectroscopic temperatures). However, the 
higher order statistics are strongly influenced by rare, 
outlying objects (e.g. major mergers), or the impact 
of observational selection effects such as cluster-cluster 
confusion. For this reason, it is unlikely that it will be 
possible to place tight priors on these parameters. 
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